Pion Deuteron Scattering and Chiral Expansions 
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We discuss convergence issues as well as short distance constraints on the deuteron wave functions 
based on chiral perturbation theory relevant to pion deuteron scattering. N on- analytical terms arise 
in the multiple scattering series for the pion-deuteron scattering length limiting the accuracy of 
the calculations. This result resembles similar findings in the structure of the NN interaction. The 
effects are found not to be numerically large. The iso-scalar irN scattering length from the iso-vector 
one and the nd scattering length yields values compatible with the experimental number and with 
much smaller errors. 
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I. INTRODUCTION 

The relevance of pion dynamics in low energy hadronic 
and nuclear processes can hardly be exaggerated 1] . Low 
energy theorems based on chiral symmetry provide a 
quantitative and model independent insight in low en- 
ergy reactions involving pions and photons and nuclc- 
ons due to the clear scale separation. The bound state 
character of nuclei makes such a description more com- 
plex from a theoretical viewpoint because many relevant 
scales enter into the problem, but one expects simplifi- 
cations to occur in the limit of weak binding and low 
energies as it is the case for the deuteron. Actually, the 
possibility of computing Pion-Deuteron scattering (for a 
review see e.g. |2J and references therein) in a model in- 
dependent way was one of the original motivations to 
introduce Effective Field Theory (EFT) approaches Q 
based on the chiral symmetry of QCD and to derive the 
corresponding low energy theorem |^. This approach 
has been vigorously extended in the last decade to the 
calculation of low energy reactions for finite nuclei (for 
comprehensive reviews see e.g. Ref. |MSI3)- This is in 
fact a rather complicated process since there are many 
corrections and physical effects which add up to the final 
result @, and it is a challenge for chiral approaches 
to make reliably calculations based on a priori estimates 
of the accuracy. The use of the original power counting 
suggested by Weinberg to extract the tt— neutron scat- 
tering length from the known 7r— proton and n— deuteron 
ones was suggested in Ref. This power count- 

ing has been modified in Ref. 11] to account for the 
long distance enhancement since the deuteron wavefunc- 
tions extend far beyond the pion Compton wavelength. 
Their result resembles the single and double scatter- 
ing terms of the widely used older approaches 0, H, 
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where the deuteron inverse moments, (r~ n ), play an es- 
sential role. Re-scattering, deuteron recoil and binding 
effects can be re-summed to all orders in a rather elegant 
formula [T^ in the isospin limit (for a review see e.g. 
Ref. 01)- Nucleon recoil effects have also been shown 
to be small [l4|. Pionic deuterium has been discussed 
in Ref. |l5j |. Traditionally, these matrix elements have 
been evaluated using potential model calculations for the 
deuteron wave functions, which generally produce diver- 
gent results for (r~ n ) (with n > 2). Perturbative wave 
functions based on OPE potentials yield divergences al- 
ready for (r" 1 ) O E3 and only recently has it been 
realized that non-perturbative calculations have a quite 
different convergent behaviour [l7L ItsL IT^ | . 

Following insightful previous works I20I l2p , w e have 
developed in a series of papers 0, l22l l23t IziL |2j| a 

framework for the treatment of the NN interaction in 
a model independent way based on long distance corre- 
lations among physical observables. In our approach the 
long distance chiral NN- One Pion Exchange (OPE) and 
Two Pion Exchange (TPE) po tentials computed within 
perturbation theory |26l 1271 l2q . are iterated to all orders 
in the Schrodinger equation very much in the spirit of the 
original Weinberg approach H although some additional 
subtleties are encountered [^jjElEil- -Due to the power 
like singular character of the chiral potentials at the ori- 
gin, a radial short distance cut-off is imposed which ul- 
timately is removed by renormalization fixing some low 
energy observables or deuteron properties. This proce- 
dure turns out not to correspond to the imposition of 
a strict power counting in the sense that corrections to 
physical observables might be estimated a priori by di- 
mensional arguments, and hence departs from the origi- 
nal EFT program. Instead, one can show that naive per- 
turbative expansions of renormalized scattering ampli- 
tudes do indeed display divergences which can be traced 
to given non-analyticities in the expansion parameter of 
the solution with the long distance potential |24| . More 
specifically, if the long distance potential is written as 
V{r) — Vope(^)+AVtpe(^)+- ■ ■ one has that the scatter- 
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ing amplitude at low energies fulfills A — -4ope = 0(X a ) 
with < a < 1 for singular potentials which fulfill the 
condition that VopeM 3> VtpeM at large distances. 
The value of the exponent a depends on the particu- 
lar structure of the potentials involved. This weird frac- 
tional power counting is uncomfortable from a theoretical 
viewpoint and admittedly non conventional but does not 
contradict the physical requirement that effects which 
are small at large distances become parametrically small 
at low energies. This happens to be so even if at short 
distances both singular potentials fulfill the opposite re- 
lation VopeC 7 ") ^ Vtpe(t). As a consequence standard 
perturbation theory cannot be applied. The singularity 
of the potentials just provides a non-analytical enhance- 
ment of the perturbation. Thus, there is a lack of system- 
atics a priori but one can always check the corrections to 
be small numerically as it turns out to be the case for the 
OPE and TPE potentials. For long distance potentials, 
which are local, as it is the case for the OPE and TPE 
potentials [26|, |23, > severe restrictions on the possible 
counterterms are derived from the physical requirement 
that the wave function be small in the short distance and 
unknown region as the cut-off is removed [24| ■ Calcula- 
tions based on finite cut-offs do not show these problems 
but then the model independence is manifestly lost since 
the particular regularization acts itself as a short distance 
model. 

With these new insights in mind we want to address in 
this paper the issue of 7rd scattering and, more specif- 
ically, the calculation of the corresponding scattering 
length within the multiple scattering formalism to all or- 
ders as sugg ested by Deloff [HI where, as already men- 
tioned, the inverse moments of the deuteron wave func- 
tion play a role. In our OPE analysis of the deuteron \Tfy 
we showed that after renormalization the first inverse mo- 
ment is finite due to the peculiar short distance behaviour 
of the deuteron wave function and because the potential 
at short distances presents a 1/r 3 singularity. The an- 
alytical structure at short distances was determined to 
high orders. A similar observation was made by Nogga 
and Hanhart Ref. |18l using the momentum space wave 
functions of Ref. |29l|. Recently, Platter and Phillips 
have analyzed this matrix element as well as the second 
inverse moment and have shown that a direct treatment 
in coordinate space allows naturally for a rather clean 
extrapolation of the renormalizcd matrix element. 

In the present work we discuss and extend these de- 
velopments to the TPE case on the light of the multiple 
scattering formalism to all orders (Sect. |nj. Let us re- 
mind that the truncation of such an expansion where the 
coefficients are in fact the inverse deuteron moments is 
motivated by the smallness of the irN scattering lengths. 
We show that higher order inverse moments are indeed 
convergent when the deuteron wave functions are built by 
fully iterating the long distance potentials (See Sect. lIII|l . 
The finite values are computed in Sect. llVl However, to 
any given approximation of the long distance potential 
there is always a finite value of the order of the moment 



above which none of them converges. In contrast, the re- 
summation of re-scattering effects provides convergent re- 
sults and indeed exhibit a logarithmic enhancement of the 
multiple scattering expansion parameter (Sect. 01. We 
use these approaches to examine the determination of the 
iso-scalar ttN scattering length within several schemes. 
Boost and finite range corrections are discussed qualita- 
tively in Sect. IVII Finally, in Sect. IVHI we come to the 
conclusions. 



II. THE MULTIPLE SCATTERING EXPANSION 
AND SHORT DISTANCE SENSITIVITY 

In a remarkable paper Deloff has analyzed ird scat- 
tering at threshold by solving the Faddeev equation and 
has worked out the scattering length within a zero range 
model based on a boundary condition for the relative 
pion-nucleon wave function. In field operator language 
this corresponds to a contact NNirir interaction, which 
can be regulated in this context by a radial regulator. In 
this approximation, intrinsic finite size effects of the ttN 
interaction are not included (see e.g. H, @ for a recent 
discussion on these and other effects). 

The result he found in Ref. [12j for the scattering 
length is quite simple 



1 + m/2M J 



dr[u{rf + w{rf] A^ d (r) (1) 



where u(r) and w(r) are S- and D-wave deuteron wave 
functions satisfying the coupled channel 3 Si — 3 Di set of 
equations 

~u"{r) + U s (r)u{r) + U sd {r)w{r) = - 7 2 w(r) , 



-w"{r)+U sd (r)u(r) 



U d (r) 



w(r) 



(2) 



and with asymptotic conditions 
u(r) -> Age-^, 
w(r) — > Ape p 



1 + 1- 



7r (7f) : 



(3) 



where 7 = V MB — 0.231605 fm 1 is the deuteron wave 
number, As is the normalization factor such that 



dr [u(r) 2 + w(rY 



(4) 



and the asymptotic D/S ratio parameter is defined by 
rj = A^/Ag. For conventions and numerical values of 
parameters we use Ref. 0, 0] throughout the paper. 

The function A^dir) in Eq. which will be called 
Deloff function for short, is given by 



~b + {b + b 1 ){b Q -2b 1 )/r 
1 -S 1 /r-(S + 6i)(6o -260/r 2 



(5) 
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with hi = (1 + m/M)bi, being 60 and b\ the irN scat- 
tering lengths according to the standard decomposition 
(assuming isospin symmetry) 



= b + bit ■ f. 



(6) 



Recoil and deuteron binding effects are taken into ac- 
count by making the simple replacements |l2^ 



h = 



(1 +m/M) 
1 + nbi 



e~ Kr /r. 



(7) 



where k = jy/2m/(m + M) = 0.117261 far 1 . The 
previous Eq. sums up all multiple scattering effects 
due to zero range nN interactions. It does not include 
Fermi motion, higher partial waves contributions nor fi- 
nite range irN corrections. We will comment on these 
corrections at the end of this paper. When expanded for 
small 60 and bi one gets the result 



0-rrd 



1 + — 

L > 2M 



6 + (6g - 26?)(J 

(bl - 2b\b - 2b\) 
(bf, - Ab\bl + 2b\) 



(8) 



which has been used quite often truncated to second or- 
der [lj. Here, the deuteron wave function average is de- 
fined as 



dr 



u(r) + w{rY 



(9) 



The multiple scattering expansion is motivated by the 
smallness of the s— wave ttN scattering lengths. Actually, 
such an expansion finds theoretical support from Chiral 
Perturbation Theory since at lowest order one has the 
Weinberg- Tomozawa (WT) relations p| 



h w 
°o 



3 [ai 
I (as 



2a 3 ) = 



Ol) 



8tt(1 + m*/M)f* 



(10) 
(11) 



which compare rather well with the exp erimental num- 
bers extracted from pionic hydrogen [3J| 



b = -(0.22 ±0.43) x 10 
&i = -(9.05 ±0.42) x 10~ 2 



(12) 
(13) 



Higher order corrections to the current algebra relations 
have been computed via standard ChPT methods |3lJ. A 
picture of the Deloff function can be seen in Fig.^for the 
WT values as well as the previous values from Ref. [30| . 
On the other hand, the measured ir~d scattering length 

is mm 



[-252 ± 5(stat.) ± 5(syst.)] x 10" 



(14) 
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FIG. 1: The Deloff function A n d(r) (in fm) as a function of 
distance in (fm). We use the Weinberg- Tomozawa values as 
well as the experimental ones for the s-wave nN scattering 
lengths and include recoil and binding corrections. 



From the viewpoint of the Chiral nuclear approach, the 
issue that we address here is to determine the compat- 
ibility of all scattering lengths using both the multiple 
scattering expansion as well as the deuteron wave func- 
tions based on chiral potentials. 

In the case of the zero range ttN interaction the con- 
vergence of the multiple scattering expansion, Eq. JHJ, 
would require in particular that any of the inverse mo- 
ments of the deuteron wave function must be finite at the 
origin. However, potential model wave functions based 
on regular potentials, i.e. r 2 U(r) — > for r — > 0, are 
dominated by the centrifugal term at short distances and 
hence satisfy the regularity conditions at the origin, 



u(r) 



w(r) 



(15) 



reflecting their L = and L = 2 angular momentum 
character respectively. So, it is clear that negative mo- 
ments fail to converge starting at third order where a log- 
arithmic divergence takes place. In addition, the short 
distance behaviour of the deuteron wave function be- 



comes relevant for 



and, in fact, potential models 



indeed exhibit this sensitivity. The third and higher in- 
verse moments are divergent. This strong short distance 
dependence looks very weird and counter intuitive, since 
we are looking at pion-deuteron scattering at zero en- 
ergy, where the wavelength of the incoming and outgoing 
pion is much larger than any of the other length scales 
of the problem. So, we regard this effect as a mathemat- 
ical artifact of the expansion and the potential model 
wave functions, and not as a genuine physical feature. 
It certainly does not agree with the philosophy under- 
lying Effective Field Theories, namely that low energy 
physics does not depend on short distance details. On 
the other hand the full formula does not present this 
problem because in spite of going to a finite limit at long 
distances, A^if) — > bo, there is a linear short distance 
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suppression of the operator, A^d{r) ~ — r, as a result of 
the re-summation and in agreement with the EFT ex- 
pectations. These considerations suggest that there may 
be problems with the convergence of multiple scattering 
expansion, which is ultimately motivated by the weak s- 
wave ttN interaction at threshold and which fits quite 
naturally within Chiral Perturbation Theory. We will 
show below that the problem is related to the implicit 
assumption of analyticity in the nN scattering lengths. 



III. SHORT DISTANCE CONSTRAINTS ON 
MATRIX ELEMENTS AND SINGULAR 
POTENTIALS 

One relevant question is whether the multiple scatter- 
ing expansion so widely used can still be undertaken and 
to what order in the case of the zero range ttN interac- 
tion. Obviously, for this to happen at a given finite order, 
say (r~ k ) < oo, we must have 



u(r) 



„{k-l)/2 + e 



, w{r) 



„(fc-l)/2 + e 



(16) 



where e > 0, and which implicitly requires that short dis- 
tances cannot be dominated by the centrifugal term. The 
only way how this may happen is that the potential be- 
comes more singular than the centrifugal barrier, which 
behaves as 1/r 2 . Indeed for a potential which diverges 
like a given power U (r) ~ r~ 2n the wave function has 
the power behaviour u(r) ~ r n l 2 (up to some exponen- 
tial or oscillatory function depending on the attractive or 
repulsive character of the potential at short distances) 1 . 
Thus, the maximal value for which the inverse moments 
converge fulfill 



< oo 



if 



U(r) 



-2k-2 + e 



(17) 



It is remarkable that chiral potentials do indeed exhibit 
these short distance singularities required by finiteness 
on the inverse moments of the deuteron wave function. 
Actually, a chiral expansion of the potential reads j2|| 

mm, 

U(r) = ^F(°)(rnr) + ^F( 2 >(mr) 



One can rewrite the expansion as 

J_,G(3) (mr)+ ... 



(18) 



(19) 



where the functions (ror) have a finite limit for van- 
ishing argument. Thus in the short distance limit 



U(r) 



M 

1 

7^6 



G<°>(0) 
G< 3 >(0) 



M 



G* 2 >(0) 



(20) 



Note that in this limit the pion mass dependence disap- 
pears 2 . In the deuteron case the coefficients in the former 
expressions become matrices [24| (see Appendix Ia")i . 

Thus, in the absence of the long distance potential none 
of the inverse moments is finite, while at LO the moments 
(r- 1 ) and (r~ 2 ) are finite, and at NLO and NNLO (r^ 1 ), 



and 



are also finite. 



Given the fact that the short distance behaviour of 
the chiral potentials does not depend on the pion mass, 
the short distance contribution to the inverse moments 
is dominated by the corresponding short distance scale, 
R, so one has 



R~ 



short 



(21) 



provided the integral is convergent. Thus, if we use the 
OPE exchange potential, R ~ M / f 2 we see that inverse 
moments do indeed become large in the limit R — > 0. As 
shown in our previous work on the deuteron [17| . long 
distance perturbation theory mistreats the behaviour of 
the wave function at short distances, introducing a very 
unnatural strong short distance dependence. Actually 
the first order contribution to the deuteron wave func- 
tion diverges. As a consequence, in the OPE potential, 
(r _1 ) also diverges as first noted in Ref. 0] using the 
PDS subtraction scheme. This was the reason to choose 
always the regular solutions of the fully iterated poten- 
tial at the origin for the OPE case 01 • The divergence 
persists also when the TPE exchange potential is treated 
in perturbation theory on the OPE distorted wave basis 
as a zeroth order approximation, since as pointed out in 
Ref. p^. the deuteron perturbative wave functions di- 
verges strongly at the origin. 

As we see, the multiple scattering expansion of the 
pion-deuteron scattering length involves negative mo- 
ments which become more convergent and hence more 
insensitive to short distance details when the NN po- 
tential is more accurately described at shorter distances. 
For this finiteness of negative moments to occur it is 
essential that the singular potentials be fully iterated 
since the short distance behaviour is highly non pertur- 
bative. This gives us some confidence on the treatment of 
the singular potentials and the renormalization process 
adopted in our previous works [Tj], |24( . As we have also 
stressed in the introduction, the successive improvements 



2 The short distance however does not coincide with the chiral 
1 In the limit n — > oo this includes also the possibility of an in- limit; in the latter case less singular subleading powers are ob- 

finiteley repulsive hard core potential. tained. 
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on the potential are parametrically small although in a 
non-analytical way for the low energy NN observables. 
Note that this is not the case for the inverse moments. 
Nevertheless, we will see that a similar non-analytical 
behaviour occurs in the perturbative treatment of the 
pion-deuteron scattering length via a multiple scattering 
expansion. 



IV. INVERSE MOMENTS FOR FULLY 
ITERATED CHIRAL POTENTIALS 

The calculation of the inverse moments for the 
deuteron wave functions for the OPE and TPE poten- 
tials is in principle straightforward. In the OPE case 
the first inverse moment was calculated in Ref. |17| and 
then in Ref. These numbers have been checked in 

Ref. I n the later reference also the second inverse 

moment (r~ 2 ) has been estimated. Here we confirm and 
extend these results to the first three finite inverse mo- 
ments in the TPE case. 

One technical aspect to consider in the present calcula- 
tion corresponds to the short distance contribution of the 
matrix elements. As an illustration we plot in Fig. H the 
integrand corresponding to the inverse square moment, 
(u 2 + w 2 )/r 2 for both the OPE and the TPE potentials. 
As we see, a substantial contribution to the integral is 
dominated by the short distance region making the con- 
vergence at short distances numerically unreliable. In 
the OPE case this can be fixed by using the analytical 
solutions found in our previous work [j"7| and computing 
the integral analytically in the short distance region. In 
the appendix lAl we analyze the problem for the OPE and 
TPE in more detail. Actually, the possibility of making 
these estimates analytically as well as the determination 
of the convergence of matrix elements is a virtue of the 
coordinate space method which exploits the locality of 
the chiral potentials in a natural way 3 . One important 
feature is that 20% of the total value comes from the re- 
gion below 0.2fm. Up to tiny oscillations the convergence 
behaves as ^/rZ for r c — > 0. 

Following Ref. we use the superposition principle 
of boundary conditions and write the deuteron wave func- 
tions as 



u(r) = u s (r) +rju D (r) , 
w(r) — ws{r) + rjwoir) . 



(22) 



For the OPE potential the short distance bevaviour has 
been displayed in Ref. ^tJ- The regularity condition at 
the origin fixes ?7opk = 0.02633. In the TPE potential 
we refer to Ref. [24( for the short distance behaviour. In 



3 This is in contrast to momentum space methods where computer 
space limitations may suggest a seeming convergence on the mo- 
mentum space cut-off but still far from the infinite cut-off limit 
result (see e.g. Fig. 5 in Ref. llSl V 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 
r(fm) 



FIG. 2: The integrand (u 2 + w 2 )/r 2 for the pion-less, OPE 
and TPE deuteron wave functions. 



such a case 77 is a free parameter which can be fixed to the 
experimental value with the corresponding uncertainties 
?7oxp = 0.0256(4). Using Eqs. (|22J) we can display the 77 
parameter dependence explicitly in observables involving 
deuteron wave functions, 



ss 



SD 



DD 



TPE 



f 



ss - 



2^1sd + ri 2 l DD 



(23) 



Numerically, for the chiral constants deduced in Ref. |3j 
(Set IV in our work) we find the expressions, 



2.43214 - 181.998 77 + 3942. rf 







\ r 1 


^TPE 


1 v 




r 2 1 


^TPE 


1 > 




r 3 / 


^TPE 



3.39582 - 189.582 77 + 4175.2 if 



3.39582- 


189.582 


v- 


- 4175.2 j] 2 


2.52146- 


202.658 




- 4350.4 ?7 2 


3.39582 - 


189.582 




- 4175.2 ?7 2 


3.19633- 


261.897 


v- 


- 5575.01 rj 2 



(24) 



And also 



e -KT 




3.08822- 


221.42577- 


f 4807.6 ?7 2 


r 


^TPE 


3.39582- 


189.582^- 


f 4175.2 ?7 2 






3.66165- 


287.55 77 + 


6189.54 ?7 2 


r 2 


^TPE 


3.39582- 


189.582 77- 


f 4175.2 ?7 2 


— 3kt 




5.14335- 


417.313/7- 


f 8914.07 t? 2 


r 3 


^TPE 


3.39582- 


189.582 77 


+ 4175.2 ?7 2 



(25) 

The numerical coefficients depend solely on the TPE po- 
tential, 7 and k for which we take the values 0.231605 and 
0.117261fm~ 1 respectively. Our results for the computed 
moments are summarized in tabic [I] Errors on the short 
distance cut-off scale linearly, quadratically and cubically 
for (r~ 3 ), (r~ 2 ) and (r^ 1 ) respectively (see appendix Ia"|| . 
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7 (fm- 1 ) 


V 






(^>(fm- 3 ) 


<^> (fm- 1 ) 


( e -^-) (fm- 2 ) 


( e -^~) (fm" 3 ) 


u(r) — e 1T 


Input 





00 


oo 


oo 


oo 


oo 


oo 


OPE 


Input 


0.02633 


0.4782(7) 


0.4208(10) 


oo 


0.380(1) 


0.3313(6) 


oo 


OPE* 


Input 


0.02555 


0.4861(10) 


0.434(3) 


oo 


0.387(1) 


0.343(4) 


oo 


TPE-Set I 


Input 


Input 


0.424(3) 


0.248(3) 


0.213(5) 


0.326(3) 


0.170(4) 


0.146(5) 


TPE-Set II 


Input 


Input 


0.424(4) 


0.248(4) 


0.214(6) 


0.327(3) 


0.171(4) 


0.147(5) 


TPE-Set III 


Input 


Input 


0.436(3) 


0.265(6) 


0.239(10) 


0.338(3) 


0.185(5) 


0.167(9) 


TPE-Set IV 


Input 


Input 


0.447(5) 


0.284(8) 


0.276(13) 


0.349(5) 


0.202(7) 


0.198(12) 


Nijmll 


0.231605 


0.02521 


0.4502 


0.2868 


oo 


0.3519 


0.2032 


oo 


Reid93 


0.231605 


0.02514 


0.4515 


0.2924 


oo 


0.3531 


0.2084 


oo 


Exp 


0.231605 


0.0256(4) 















TABLE I: Inverse moments of the radius and e~ Kr jr moments for the deuteron. We consider the OPE and TPE potentials; in 
the case of the OPE potential we have taken g-nNN = 13.08 (OPE) and qa = 1.26(OPE*), while in the TPE case we show the 
results corresponding to the four set of chiral couplings considered along this work. In the OPE case the error is estimated by 
varying the semiclassical matching radius in the range 0.1 — 0.2 fm (see Appendix J^), while in the TPE case the error comes 
from the experimental uncertainty of the D/S ratio, r\ — 0.0256(4). TPE Sets I,II,II and IV refer to the chiral parameters, ci, 
C3 and C4 of Refs. [^], [23],[2i| and [3] respectively. 



Besides the fact that the non-perturbative inclusion of 
long distance chiral potentials to all orders provides con- 
vergent inverse moments, one of the features one observes 
from inspection of table [I] is the size of the changes in- 
duced when going from OPE to TPE potentials. Whereas 
the TPE effects have moderate effects on the (r -1 ) mo- 
ment as compared to the OPE value, almost a factor 
of two reduction for (r~ 2 ) is obtained. This is to be 
expected as the TPE potential modifies the short dis- 
tance region 4 . Moreover, we get a large reduction of 
the second inverse moment since (?"- 2 )ope = 0.41fm -2 



while 



}tpe = 0.25 — 0.28fm very close to the value 



0.286 - 0.345fnr quoted in Ref. [nLAs we also found 
in our previous work [24| the Set IV |34j] agrees best with 
the Nijmll and Rcid93 potential values. 



V. DIVERGENT INVERSE MOMENTS AND 
NON-ANALYTICAL BEHAVIOUR 

As we have also discussed, the use of the Deloff func- 
tion (with or without binding and recoil) produces finite 
numbers regardless on the short distance behaviour of the 
deuteron wave function assuming it can be normalized. 
On the other hand when the multiple scattering expan- 
sion is undertaken divergences appear at some stage. Al- 
though this appears to be a bit puzzling, the mathemat- 
ical reason why this is happening is indeed quite simple. 
If we use for definiteness the WT values for the wN scat- 
tering lengths, the denominator has complex conjugated 
poles located 



(l±n/7> 
16tt/ 2 



(26) 



Finally, let us note that, as already argued above, 
the more singular the potential at short distances the 
more convergent the matrix elements, and in particu- 
lar short distance contributions to the TPE potential 
are much smaller than those in the OPE case. For in- 
stance for the second inverse moment in the region be- 
low r c = 0.2fm one gets (r- 2 )^°g = 0.08fm- 2 while 
(r- 2 )f°£ = 0.003fm- 2 . This also suggests that finite 
cut-off effects in matrix elements become smaller when 
the long distance potential is improved at lower distances. 



4 A similar factor of two was also found in the effective range 
parameter for the 1 So phase shift when going from OPE to TPE. 



where binding and recoil corrections have been neglected 
for clarity. In this approximation the radius of conver- 
gence becomes |r| = m/47r\/2/ 2 which is a rather small 
distance and goes to zero in the chiral limit. For finite 
values of b\ one can only expand above that region. The 
relevant scale is given by b\ ~ — 0.2fm so that in the limit 
of small b\ we become sensitive to the short distance be- 
haviour. Numerically the radius of convergence is given 
by \t c \ = 0.1 5fm, so the the large r expansion converges 
only for r > \r c \. Thus, if we cut-off the integrand be- 
low such a value we get a convergent expansion but at 



5 The full expression is 
\re Kr \ 



y/2m 



(1 + m/M) [8?r/ 2 (l + m/M) - mn] 
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the same time an important piece of physics is neglected. 
The short distance enhancement can be seen by display- 
ing at the function A^(r) in Fig. ^ The bump takes 
place at r ~ 0.2fm so it is indeed true that a relevant 
contribution comes from short distances. 

The kind of non-analyticity that appears can be easily 
illustrated with the asymptotic deuteron wave function 

u(r) = v / 27e" 7r , w(r) = , (27) 

which corresponds to the pion-less theory. Even though 
the first term in the multiple scattering expansion in 
Eq. |JSJ does not converge, the full formula given by 
Eq. (JTJ yields a finite analytical result. Neglecting m w /M 
corrections and taking bo = for illustration purposes 
one gets 

a* d = 4& 2 7 /(-6 l7 ) (28) 

where 



J (x- xi)(x - x 2 ) 
e- ta; ^ir(0, -fxi) - e- te ^ 2 r(o, -tx 2 ) 
x\ - x 2 

In the case bo = we have x\^ 2 — (1 ± iv7)/2 and 
F(0, z) = J e~ t /tdt is the incomplete gamma function 
which has a logarithmic branch cut at z = 0. Thus, there 
is a branch cut singularity at b\ = 0, the leading term 
being 

a7ld 4^ 7 log(-6i7) = -3.06 x lCT 2 ™^ 1 (30) 

This example shows explicitly the kind of enhancement 
that one might expect. The full result with the short 
distance wave function and using the WT nN scattering 
lengths a^d = —2.815 • IO^to^ 1 , a quite reasonable value 
taking into account the poor quality of the wave function. 
This value overshoots the real value due to the fact that 



the short range wave function, Eq. I|27(l does not vanish 
at the origin. For more regular functions such as those 
of potential models we expect that logarithmic short dis- 
tance enhancement takes place precisely at the fourth 
and higher orders where the divergence of the inverse 
moment becomes manifest. The chiral TPE deuteron 
wave functions present the logarithmic enhancement at 
fifth order. 

In Fig. |3| we plot the integrand for a variety of wave 
functions. A good feature of the use of the Dcloff func- 
tion is the irrelevance of short distance behaviour as com- 
pared to the multiple scattering expansion where there is 
an enhancement of short distances. One important les- 
son we learn from this exercise is that the appearance of 
non-analytical enhancements found in non-perturbative 
treatments of the NN interaction extends also to the com- 
putation of matrix elements. Strict power counting sim- 
ply does not hold, although the corrections are paramet- 
rically small. Numerical results for the pion-less, OPE 
and TPE potential cases using the Dcloff formula are 
presented in Table [H] for the central values of the ttN 
scattering lengths. 

As we see, the multiple scattering series provides can- 
cellations which are independent on the short distance 
constraints for the inverse moments discussed in Sect. lIIll 
One might think that they might be correlated through 
the chiral expansion, in the sense of a perturbative re- 
ordering of the truncated multiple scattering series by 
re-expanding the convergent inverse moments. This is 
unlikely, since the kind of non-analyticities appearing in 
the multiple scattering series as a function of bo and bi 
and those in the inverse moments as a function of the 
chiral potential parameters are of quite different nature. 

It is worth displaying numerically the convergence of 
the multiple scattering scries for both the OPE and the 
TPE deuteron wave functions. For illustration purposes 
we have taken the Dcloff function ignoring both binding 
and recoil corrections (similar features are observed if 
these corrections are taken into account). We obtain 



a^lopE = 130.42 6o + 62.32 (^-2 6?) + 54.84(6g-2 6 6i-26j)+0(6 4 log6) (31) 

-0.47013 -2.744 0.367 

q^Itpe = 130.42 bp + 58.25 (bj - 2 1\) + 37.01 (pi - 2 b b\ - 2b\) + 35.97 (feg - 4 b%b\ + 2bj) +Q{b 5 log b) , (32) 

-0.47013 -2.565 0.248 0.035 



where the results given in the under-braces are in units of 
lO -2 ™^ 1 (as the pion-deuteron scattering length), and 
the biS are in fm. As we see, the bulk of the contribution 
is given by the double scattering term. Using the full 
OPE and TPE potentials and the full A^dir) operator 



with no recoil and binding we get 

a° d PE = -2.873(3) • IO^to" 1 (33) 
a TPE = -2.77(2). l0- 2 m-\ (34) 

while the analytic contributions in Eq. 132|) sum up to 
—2.847 and —2.75 respectively. As we see, the non- 
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FIG. 3: The integrand {u 2 + w 2 )A^ d {r) for the pion-less, OPE 
and TPE deuteron wave functions. 

analytical effects are not dramatic, which was not com- 
pletely obvious a priori and to a certain extend are com- 
patible with the uncertainties in the TPE case. Obvi- 
ously, very precise estimates might be sensitive to the 
non-analytical pieces. In any case, it would be rather 
interesting to compute the non-analytical contributions 
per se. 

Finally, we plot in Fig. 21 the correlation between the 
iso-scalar and iso-vector irN scattering lengths b and b\ 
deduced from direct application of the Deloff formula in 
several schemes where errors from the n — d scattering 
length and chiral potential parameters are taken into ac- 
count. The main source of uncertainty in all cases turns 
out to be a^d- If we take the iso-vector b\ value from 
Ref. |3(j we obtain in the TPE case including both bind- 
ing and recoil corrections the result 

b Q = -0.3(1) • lO^m" 1 , (35) 

which is compatible with the measured value but about 
an order of magnitude more accurate. 



VI. REMARKS ON BOOST AND FINITE 
RANGE CORRECTIONS 

One of the effects we have not taken into account in 
our discussions has to do with the boost corrections due 
to the fact that the CM it — d and it — N systems do 
not coincide, so the nucleons inside the deuteron recoil 
differently as the deuteron. Although it has been argued 
that these effects are small 14|] it is interesting to rean- 
alyze the issue on the light of the present investigation. 
Finite range ttN corrections have also been computed in 
Refs. U and the leading contribution involves similar 
operators as in the boost corrections case, so we will refer 
mainly to the latter case in the following. Actually, the 
effect can be estimated perturbatively yielding an 0(Q 4 ) 



FIG. 4: The range of values of the nN scattering lengths &o 
and &i in units of lO -2 ™,^ 1 that reproduce the pion-deuteron 
scattering length for different approximations (see main text). 
The error bars are generated from the experimental error in 
a-wd- We use the TPE deuteron wave functions with r\ = 
0.0256 and Set IV Hi. 





Static 


Recoil & Binding 


u(r) — e 7r 


-4.605 


-4.313 


OPE 


-2.873(3) 


-2.446(3) 


OPE* 


-2.905(3) 


-2.480(2) 


TPE-Set I 


-2.676(11) 


-2.238(12) 


TPE- Set II 


-2.675(12) 


-2.237(12) 


TPE-Set III 


-2.725(16) 


-2.29(2) 


TPE-Set IV 


-2.77(2) 


-2.34(2) 


Nijmll 


-2.786 


-2.354 


Reid93 


-2.789 


-2.357 



TABLE II: Pion-deuteron scattering length in units of 
10 -2 m" 1 for OPE and TPE potentials in different approxi- 
mations. Static means the Deloff formula without recoil or 
binding, and recoil & binding taking into account the re- 
coil and binding corrections according to Eq. J7J. We take 
b = -0.22 • 10 -2 m" 1 and 6i = -9.05 ■ 10~ 2 m' 1 . OPE and 
OPE* account for taking q^nn = 13.08 and qa = 1-26 re- 
spectively as input. In the OPE case the error is estimated 
by varying the semiclassical matching radius in the range 
0.1 — 0.2 fm, while in the TPE case the error comes from the 
experimental uncertainty of the D/S ratio, r\ = 0.0256(4). 
TPE Sets I,II,II and IV refer to the chiral parameters, ci, C3 
and c 4 of Refs. |2j|, B-EIJl and respectively. 



correction according to the modified counting proposed 
in Ref. HJ, 



|Boost = - l + m,/2M87rM3/, a (gi " 8MC2){P2) 

(36) 
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where 



(P 2 ) = 



dr 



u'(r) 2 +w'{rf +6 



(37) 



The operator (p 2 ) is weakly non-local and strictly speak- 
ing this number is infinite for any singular potential. Ac- 
tually, the divergence behaves as l/y/r^ for the OPE po- 
tential and as 1/rj! for the TPE potential. The large 
values of (p 2 ) have been discussed in finite cut-off cal- 
culations in momentum space and attributed to the in- 
ner knots of the wave function, hence favoring the use 
of relatively small cut-offs. This appears to be a serious 
difficulty, but resembles the one faced already for the in- 
verse moments; one expects non-analytical effects in the 
couplings in front of (p 2 ). This suggests further lines 
of research, in particular invoking physically motivated 
boost re-summations or including nonlocal effects in the 
deuteron wave functions. As it is known, the higher or- 
der contributions to the NN potential, 0(Q 4 ) according 
to Weinberg's counting, contain nonlocal pieces. In mo- 
mentum space and up to NNLO the long distance part 
of the potential depends on the momentum transfer q 
only and not on the total momentum k. Essential non- 
localities, i.e. contributions of the form V(q, k) = L(q)k 2 
with L(q) a non-polynomial function, depend weakly on 
the total momentum and appear first at N LO [36| i.e. 
also 0{Q A ) due to relativistic 1/M 2 one loop contribu- 
tions. In coordinate space this weak non- locality corre- 
sponds to a modification of the kinetic energy term in the 
form of a general self-adjoint Sturm-Liouville operator, 
—u"(r) — ► — (p(r)u' (r))' , with a singular p(r) function at 
the origin and exponentially decaying at long distances. 
It is at present unclear how these non-locally modified 
wave functions might influence the pion-deuteron scat- 
tering length and further work along these lines should 
be pursued. A more promising perspective consists of 
resumming boost corrections non-perturbatively 6 . Ob- 
viously, to achieve a definite conclusion on this issue one 
should reexamine the nNN system within a full quantum 
mechanical Faddeev approach including chiral potentials 
between the nucleons which is nontrivial. We leave such 
a study for future developments. 



has the important feature of providing a physically mo- 
tivated suppression of the deuteron matrix elements at 
the origin. This is most welcome since it agrees quali- 
tatively with the naive expectation that low energy pro- 
cesses should not depend strongly on short distance de- 
tails. The finiteness of the multiple scattering expansion 
of the full expression truncated to a given order requires 
a series of short distance constraints which are remark- 
ably fulfilled by the regular solutions of the deuteron 
wave functions corresponding to chiral potentials. This is 
so precisely because the potentials become singular and 
are iterated to all orders. Nevertheless, further insen- 
sitivity at short distance can be gained by using a re- 
summation formula proposed by Deloff. We have noted 
that divergences in multiple scattering expansion arise 
because non-analyticities in the chiral expansion appear. 
This is a rather general feature which goes beyond just 
the example of pion-deuteron scattering addressed in this 
paper and will also occur in other low energy reactions. 
The found non-analyticitics are not large numerically due 
to the good short distance behaviour of the deuteron 
wave functions. Finally, we have compared several ways 
of extracting the poorly known iso-scalar ttN scatter- 
ing lengths from the known iso-vector ttN and the pion- 
deuteron scattering length, yielding, as expected, a com- 
patible but more accurate result. 
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VII. CONCLUSIONS 

In this paper we have analyzed a particular re- 
summation of re-scattering effects to pion-deuteron scat- 
tering length corresponding to the fixed center approx- 
imation and where recoil and binding effects may be 
taken into account to all orders. This re-summation 



For instance, for the s-wave ttN term one would have the replace- 
ment 0: mt(p 2 } -* (p 2 ( A f+ZrV-p^ )- This corres P onds to a 
finite nucleon recoil, which has better converging properties. 



APPENDIX A: SHORT DISTANCE 
CONTRIBUTIONS OF MATRIX ELEMENTS 
INVOLVING INVERSE MOMENTS 



In this appendix we quote analytical results for the 
short distance contribution to the inverse matrix ele- 
ments. Most integrals are sufficiently converging by let- 
ting the short distance cut-off to go to zero. The 1/r 2 for 
the OPE case is a bit special since the integrand presents 
many oscillations which make numerical extrapolation 
unreliable. As pointed out in Ref. [l^ it is much bet- 
ter to use the analytical expressions deduced in Ref. ^} • 
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Short distance OPE 



where 



At short distances the OPE potential behaves as 



U? PE (r) - ^ 



U? d PE (r) 
U OPE 



(r) 



R s d 

Rd 



where R d = AR, R sd = 2y/2R and R = 3g 2 A M/32irf 2 = 
1.07764fm. Going to the diagonal basis the solution can 
be written as 



u(r) 



^v A (r) - -j=v R (r), 



w(r) -► -^=v A (r) + sj ^v R (r) , 



(A2) 



Mr) = 



>' \ 



3/4 



-4VV? 



2i?e 



(A3) 



(Al) v A {r) = ( 



r v 3/4 
£7 



The constants Ci#, C 2 fl, Cia and C 2 a have been fixed 
in Ref. [ljj from matching the numerical solution to a 
short distance expansion (the regularity condition Cm = 
is imposed). The integral can be computed for short 
distances yielding 



R 2 [ r ° , u 2 +w 2 



OPE 



384^1 



W2 ClA ° 2A { 
384 ° 2A i 



91264 tt+ 1536^+ 140 x| + 22816 y/x~ c cos(——) 



1085x! cos(^L) + 2660x c sin(-^L) + 182528 Si(—^=) } 

\ / Xr \/ 30 r \/Xr ' 



2660 x c cos(-^L) + 182528 Ci(-)L) + 31 (-736 + 35 x c ) sin(-)L)} 



91264 tt + 1536 ^/x~ c + 140 x c 2 - 22816 ^/x~ c cos( — ) 

\l Xr 



1085x1 cos(^L) - 2660 x c sin(-^L) - 182528 Si(-)L))} + 0(x 2 c ) (A4) 

Xq \/Xr- \/Xr- * 



with x c — r c /R and Si and Ci are the sine and cosine behaves as as |26l 1271 126 
integral functions respectively defined as 



Si(z) = 
&(z) = 



sint 



dt 



sint 



-dt 



(A5) 



The convergence on r c is shown in Fig. |SJ We see that wnere 
0.06fm" 2 , i.e. about 20% of the result, comes from the 
region below 0.2fm. We have checked that up to this 
region the result is rather stable with the given terms. 



Short distance TPE 



(Rs) 4 
(R sd ) 4 
(Rd) 4 



Uj PE (r) 
Uj d PE (r) 
Uj PE (r) 



3.9 



Rj 

R td 

Rd 

r 6 



(A6) 



128/> 
3V2g 2 



^(4-3 3 i + 24c3-8c 4 ) 



4 (-4 + 3g 2 A - 16c 4 ) 



128/> 2 

(-l + 2g 2 A + 2c 3 -2c 4 ) (A7) 



9g 2 

32/47T 2 



The short distance behaviour of the TPE has been de- and a = Mci are the low energy chiral couplings appear- 
termined in Ref. |24| . The potential at short distances ing in ttN scattering. Diagonalizing the corresponding 



11 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

r c (fm) 



FIG. 5: The short distance cut-off dependence (in fm -2 ) of 
the integrated negative second moment from r c to infinity, 
J^°(u 2 + w 2 )/r 2 dr as a function of r c (in fm) 



matrix 



i?„ R' 



d - 



Rt 



R d 



cos a 


sin 


— sin 


9 cos 


cos 9 


— sin 


sin 9 


cos 6 




-R 4 



(A8) 



Note that the potential is negative definite. In the diag- 
onal basis one has 




/ cos 9 sin 9 
I — sin 9 cos 

where the short distance eigen functions are 
v+{r) = 
V-(r) = 



(A9) 



r 



r 

~R~ 



C+, s sin 
C_ , sin 



lRl 
2 r 2 

\Rl_ 
2 r 2 



C 4 



C- 



cos 



-1R%- 




2 r 2 


} 


-\R 2 _- 


1 


2 r 2 





(A10) 



Higher order contributions could in principle be com- 
puted similarly to what was done in the OPE case 0] al- 
though such refinements will not be needed here. Match- 
ing our numerical solutions to these short distance solu- 
tions at r c — O.lfm we get (for Set IV of parameters) 



R + = 2.11fm and R- = 1.16fm 

C+, c = 6.174- 257.329?7 

C*+, s = -1.585 + 38.745r; 

C-. c = 2.451 - 118.8477 

C_ iS = 2.522 - 84.4077 



(All) 



for the long distance normalization such that u(r) 




can be 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

r c (fm) 



FIG. 6: The short distance cut-off dependence of (in fm -3 ) 
the integrated negative second moment from r c to infinity, 
J^°(u 2 + w 2 )/r 3 dr in the TPE case as a function of r c (in 
fm). We use Set IV of chiral constants. 



computed analytically using that 



w 



dr = 



-dr 



(A12) 



and the appropriate normalization. For the third inverse 
moments the integral can be written in terms of Fresnel 
integrals. The result is depicted in Fig. As we see, 
the short distance contribution for the moments is rather 
small. The contributions from the region between the 
origin and r c = 0.2fm are given by 



/„-i\short 

\' /TPE 
/ -2\short 
\' /TPE 
/ -3\ short 
\' /TPE 



0.0004fm 
0.003fm 
0.025fm 



-l 

-2 
3 



(A13) 
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